A lattice calculation of vector meson couplings to the vector and tensor currents using 

chirally improved fermions 
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We present a quenched lattice calculation of fy / fv, the coupling of vector mesons to the tensor 
current normalized by the vector meson decay constant. The chirally improved lattice Dirac opera- 
tor, which allows us to reach small quark masses, is used. We put emphasis on analyzing the quark 
mass dependence of fy / fy and find only a rather weak dependence. Our results at the p and <f> 
masses agree well with QCD sum rule calculations and those from previous lattice studies. 



I. INTRODUCTION 

It is widely accepted that the rich spin structure of 
hard exclusive processes involving light vector p,K*,(f> 
mesons provides a number of nontrivial possibilities to 
study the underlying short-distance dynamics. One ex- 
ample is given by vector meson electroproduction at 
large virtualities or large momentum transfers at HERA 
0,1113, HERMES 0, and, in the future, COMPASS. 
In these cases the theoretical predictions for the produc- 
tion of longitudinally and transversely polarized mesons 
are very different |5j. Other examples are provided by 
exclusive semileptonic B — > ptvg, rare radiative B — > fry 
or nonleptonic, B — > np etc. decays of S-mesons, which 
are attracting continuous interest as the prime source 
of information about the CKM mixing matrix; see, e.g., 
Ref. [|| for an exposition of recent developments. In all 
cases, disentangling the longitudinally and transversely 
polarized (final) vector-meson states proves to be crucial 
since these cases often correspond to different underlying 
weak interaction physics. The theoretical description of 
such processes is developing rapidly and thus requires a 
more accurate and reliable determination of the relevant 
nonperturbative parameters. 

From the vector meson side, the QCD description in- 
volves vector meson distribution amplitudes 0, @ 
which correspond to probability amplitudes for finding 
the quark and the antiquark in the meson with given 
momentum fractions and a small transverse separation. 
The distribution amplitudes for longitudinally polarized 
and transversely polarized vector mesons are different 
and their normalization (i.e. the integral over the mo- 
mentum fractions) is given by the matrix elements of the 
vector and the tensor current 
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respectively. Here V(p; A) is a generic light vector me- 
son with momentum p and polarization vector e? such 
that e%p a = 0, p 2 = m 2 v . Also o a(i = §[7° ,Y] and 
q is the light quark field, q = u,d,s . . . . In Eqs. 0J, 
(|2"|l we suppress the isospin structure, for brevity. The 



vector couplings /y are known from the experimental 
measurements in leptonic decays |10| . while the tensor 
couplings fy have to be calculated using some nonper- 
turbative approach. In particular, from QCD sum rules 
[TT| one obtains 



f^ = 160±10MeV, 



(3) 



(157±5 MeV in Ref. Qjg) and in the applications usually 
a very weak dependence on the quark mass of the ratios 
was assumed 03, [Dl 



f 1 - 
fp 



f K * 



(4) 



(A n earlier QCD sum rule determination of fp in Refs. 0, 
IhH ] suffers from a sign error, see Ref. .) At this place it 
is necessary to add that the couplings fy, in contrast to 
fy, are scale dependent. The corresponding anomalous 
dimension is equal to [TEl Irrl Il7| 
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(The three-loop terms [l8|, left out above simply for 
brevity, arc included in our analysis.) The number in 
© is given for a low normalization point of p — 1 GeV. 

The error given in © does not include intrinsic un- 
certainties of the sum rule method itself which are diffi- 
cult to quantify. Therefore, in view of the importance of 
this nonperturbative input for phenomenology, an inde- 
pendent confirmation of these numbers in a lattice cal- 
culation with controllable errors is extremely welcome. 
Earlier attempts were presented in Refs. [l!j|2(| and the 
results agree well with the numbers from sum rules (there 
is also a current effort by the QCDSF Collaboration to 
determine /r~; the preliminary result also agrees with the 
sum rule result 21]). These calculations were done with 
0(a) improved Wilson fermions. With this choice, how- 
ever, one is restricted to relatively heavy pseudoscalars 
and the smallest pseudoscalar-mass to vector-mass ratio 
reached in [2(j is mps/my = 0.56. 
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Here we present a new calculation of the matrix ele- 
ments which is complementary to the first lattice results 
|2fij . We use the recently developed Chirally Improved 
(CI) Dirac operator |22ll23| . It is based on the Ginsparg- 
Wilson relation |2J] and has been shown |25|, l26l 127] to 
allow simulations with pseudoscalar-mass to vector-mass 
ratios down to mps/my = 0.28 at relatively small cost. 

In the results we present here we include lattice data 
down to mps/mv = 0.33 and thus push our calcula- 
tion considerably closer to the chiral limit. Also at the 
heavy end we have additional data points such that our 
quark masses cover a range of mps/mv more than twice 
as large as the range in J2jj. This allows us to check 
the weak quark-mass dependence of fy/fv predicted by 
QCD sum rules. The data are compatible with a linear 
behavior in the quark mass, respectively in mp s , and the 
slope is small. Only for our lightest quarks we find a de- 
viation, which is, however, a finite size effect. Our final 
numbers agree well both with the QCD sum rule results 
and the lattice calculations in Refs. [L9II20I . 



II. PARAMETERS OF THE CALCULATION 

We perform a quenched calculation with configurations 
from the Liischer-Weisz gauge action [28, 29] and one step 
of HYP blocking [3(j ■ The final numbers we quote were 
computed on 16 3 x 32 lattices at two different values of the 
gauge coupling, giving rise to lattice spacings of a = 0.15 
and a = 0.10 fm [23|. For these lattices the parameters 
are listed in Tabled In addition we performed a series of 
tests on smaller lattices (12 3 x 24, 8 3 x 24) and at larger 
cutoff (a = 0.08 fm). This serves to analyze the influence 
of finite volume and to study the scaling behavior. 

We compute fermion propagators using the CI Dirac 
operator at 11 different bare quark mass values ranging 
from am = 0.01 to am = 0.20. These quark masses cover 
a range of mps/my = 0.33 to mps/my = 0.92. We do 
not encounter exceptional configurations and could work 
at even smaller quark masses. However, as we will show, 
decreasing the quark mass further, thus going closer to 
the chiral limit, is not sensible for a calculation of the 
matrix elements at the current physical volumes due to 
finite size effects. 

In [25H26II27I I it is demonstrated that the CI operator is 
effectively 0(a) improved. Furthermore the spectroscopy 
results presented there show only a very small variation 
in a 2 and thus are essentially free of scale dependence. 
However, improved operators are not yet available for the 
CI operator and thus for the matrix elements computed 
here one can expect only linear scaling in a for our final 
results. 

The renormalization constants necessary for our ob- 
servables were calculated according to the method pre- 
sented in [22, H3 • The procedure is patterned after the 
definition used in the continuum. In a fixed gauge (we 
use the Landau gauge) the numerically evaluated ampu- 
tated Green's functions are compared to their tree level 



counterparts and the renormalization constants are read 
off. The resulting numbers are in the so-called RI-MOM 
scheme and are converted to MS using the perturbative 
matching coefficients. For the vector and tensor renor- 
malization constants, however, these matching constants 
are 1 up to NLO in perturbation theory. A detailed anal- 
ysis of renormalization for the CI operator will be pre- 
sented in [3 1| and here we only quote the numbers we 
need. 

The renormalization constants were evaluated for the 
same set of quark masses that was also used in the calcu- 
lation of the matrix elements. When plotted as a function 
of the quark mass, the data are found to follow a straight 
line very well and the chiral extrapolation is straightfor- 
ward. This procedure was repeated for several values 
of the 4-momentum p. Following |20| we base our cal- 
culation on the numbers extracted at the cutoff, i.e. at 
p 2 = a~ 2 . To get the numbers at exactly a~ 2 we linearly 
interpolated the chirally extrapolated values of Zy , Zt 
between the two momenta with p 2 just above and below 
a~ 2 . The resulting numbers for Zt were then evolved to 
our target scale of [i = 2 GeV using the renormalization 
group equation. Zy remains constant as a function of 
/j, (no anomalous dimension) and we used the value at 
fi 2 = a~ 2 . Our numbers for Zy and Zt are listed in 
Table CH 

III. EXTRACTION OF THE RAW DATA 

If we only consider the vector-meson contribution to 
the vector-vector (VV) and vector-tensor (VT) correla- 
tors and contract two of the Lorentz indices, we arrive at 
the following expressions: 

9^ (0 \q^q\ V(p; A)) (V(p; A) \ffq\ 0) = 

m V fv ( e A e A*S^) . ( 6 ) 

9v» (0 197^1 V(p; A)} (V(p; A) \qa av q\ 0) = 

im v f v fy-(n)p a (e£er<&«0 - (?) 

where e l te^*g fJ , ll = —3. Taking the ratio of these (vector- 
tensor/vector- vector), we have 

R a = . (8) 

m v fy 

with one remaining Lorentz index from the tensor cur- 
rent. It is this ratio we wish to extract in our lattice 
calculation. 



p 


a (fm) [21] 


L (fm) # confs. 


am q 


m PS /m v [26J 


7.90 


0.15 


2.4 100 


0.02 - 0.20 


0.38 - 0.85 


8.35 


0.10 


1.6 100 


0.01 - 0.20 


0.33 - 0.92 



TABLE I: Parameters for our large lattices (16 3 x 32). 
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p 


Zvia- 1 ) 


Zria- 1 ) 


Z T (2 GeV) 


Z T /Z v (2GeV) 


7.90 


0.9346(7) 


1.0598(10) 


1.0247(10) 


1.0964(19) 


8.35 


0.9780(14) 


1.0542(19) 


1.0532(20) 


1.0768(37) 



TABLE II: Results for the renormalization constants. The 
numbers in parenthesis give the scale fi. It is either the cutoff 
{/! = a- 1 ) or n = 2 GeV. 



For the non-zero-momentum correlators, however, 
there is an additional complication, which arises because 
the quark sources are smeared in the spatial directions. 
Since the smearing is not performed in the time direction 
as well, Lorentz symmetry is lost, affecting the e v \p v = 
condition (checking this relation explicitly with the differ- 
ent Lorentz components of the 0, VT correlators, we 
find the discrepancy to be small, but significant). More 
explicitly, we have for the VT correlator 

(0 \(Uq)r(Uq)\ V(p; A)) (V(p; A) \qa a »q\ 0) = 

im v f v tifji) Zu(p) e» (eTp a - ef*p») , (9) 

where U represents the smearing operator and Zjj(p) is 
the corresponding factor for the amplitude. For the a = 
j case, we may no longer sum over all of the Lorentz 
indices, fj, — v, and still expect the second term to vanish. 
We may, however, consider only the contributions where 
fi = v = k 7^ j and p k = 0. Then, using only the 
v = k, p k = contributions for the VV correlator as 
well, we have the appropriate cancellations in the ratio 
R? (p) . This finite-momentum ratio yields much larger 
statistical errors than the zero- momentum ratio, R (0), 
and we pursue this approach only as a consistency check. 

We note here that the smearing has no effect upon the 
zero-momentum results since these have contributions 
only from the spatial components of the vector-current 
source and the smearing amplitude, Zu(p), is the same 
for all spatial directions at p = 0. 

The Euclidean lattice formulation provides a natural 
selection of the desired vector-meson contribution to our 
correlators; we simply need to ensure that our two cur- 
rents are separated by a large enough distance in the 
imaginary time direction. The vector-vector correlator, 
which corresponds to Eq. @ , thus reads 

C vv (p;t) = £ ((^)o,o(^0*, t ) ' ( 10 ) 

X 

We project to definite momentum p, using the necessary 
phase factor, and sum over the relevant, spatial Lorentz 
indices, k (pk = 0). For the tensor current we have to 
distinguish two different cases: When the open tensor 
index is time (= 4) we consider the temporal correlator 
C\ T defined as (k is summed) 

cx T {p\t) = ((^7fcV0 o ,o(^^W) ■ 

X 

(11) 



For open spatial tensor indices j = 1,2,3 we study the 
spatial correlator Cj T defined as 

Cy T (p;t) = £ e** ((^0 Ol o(^)*, t ) < 

X 

(12) 

where pj ^ and the sum over k is restricted to the two 
values different from j, such that pk = 0. The tempo- 
ral and spatial correlators C\ T and Cj T differ in two 
respects: Firstly, the spatial correlator is non-vanishing 
only for non-zero momenta, while the temporal correlator 
gives a contribution also at zero momentum. Secondly, 
the spatial and temporal correlators differ by their time 
reversal symmetries. In particular one finds: 

C vv (p-t) = A vv (p) (e- E vW + e Ev(m-T)^ + 

CX T {p;t) = A v J(p) (er E ^ - e EAp){t-T)^ + 

Cj T {p;t) = Aj T (p) (e- E vtv)t + e E v (v)(t~T)^ + 

(13) 

The dots indicate corrections due to excited states which 
play a role only at small time separations but become 
suppressed exponentially for larger values of t. 

Let us have a first look at our data for C vv and C VT , 
separately. This serves to check how well the expected 
behavior of Eqs. I|13|) is seen and how serious the effects 
of excited states are. For this test and also for the ex- 
traction of the data we fold all correlators about T/2 
according to their time-reversal symmetries and average 
them appropriately. In the plots we thus show only the 
folded and averaged correlator in the range t = ... T/2. 

In Fig. we show effective mass plots, i.e. 
ln[C(t)/C(t + 1)], for our folded correlators defined in 
Eqs. $TUb - (H2J). The effective mass E vv for the vector- 
vector correlator (top plot) and the effective mass E\ T 
for the temporal vector-tensor correlator (middle plot) 
are shown for zero momentum, while the effective mass 
Ej T for the spatial vector-tensor correlator (bottom 
plot) is shown for momentum p 2 = (2ir) 2 /L 2 . The 
data for the figure are from the 16 3 x 32 lattices with 
a = 0.15 fm and two values of the bare quark mass, 
am q = 0.05,0.10. The spatial correlator has been aver- 
aged over all three components j = 1,2,3. However, the 
scatter of the results for the individual components j is 
small. 

Both vector-tensor correlators show plateaus between 
t/a = 4 and t/a = 12 and the vector- vector correlator 
for t/a — 6 ... 12. For smaller time separations it is obvi- 
ous that contributions from excited states are important. 
For t/a larger than 12 we observe two effects: Firstly, an 
increase of the error bars which is simply an effect of the 
reduced signal to noise ratio for longer propagation in t 
and, secondly, a systematic deviation from the plateau 
due to the sinh / cosh behavior of the correlators. For all 
quark masses we looked at we were able to identify stable 
plateaus. When comparing our different results for the 
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FIG. 1: Effective masses for the vector-vector (top plot), tem- 
poral vector-tensor (middle plot) and spatial vector-tensor 
(bottom plot) correlators. The effective mass for the spa- 
tial vector-tensor correlator is shown for momentum p 2 — 
(2tv) 2 /L 2 , while the other two effective masses are for zero 
momentum. The data shown here are for the 16 3 x 32 lat- 
tices with a — 0.15 fm and two values of the bare quark mass, 
am q = 0.05,0.10. 



vector meson energy Ey{p) which determines the time 
behavior of all three correlators [see Eq. I|13|)] we found 
that the values of Eyip) agreed well for the different 
operators. Note that the spatial vector-tensor correlator 
shown in the bottom plot of Fig. ^ is for nonzero momen- 
tum and thus its plateau is slightly higher than for the 
other two correlators. 

Once the amplitudes are known we can form the tem- 
poral and spatial ratios R4 and Rj defined as 



R 4 (p) = A v J{p)/A vv (p) 



Rj(p) = A] T {p)/A vv {p) 



vv I 



(14) 



We have already remarked that R4 is non- vanishing also 
at zero momentum, while Rj gives a contribution only 
for non-zero momenta. Since the quality of lattice cor- 
relators generally is decreasing with increasing momenta 
we evaluate i?4 at zero momentum and Rj at the small- 
est possible momentum, i.e., p 2 — (2ir) 2 /L 2 . The ratios 



a (fm) 


am q 


t vv 


t VT 


X 2 /d.o.f. 


i? 4 (0) 


0.15 


0.02 


8-12 


4-12 


10/11 


0.678(36) 


" 


0.03 


8-13 


4-13 


11/13 


0.680(24) 


" 


0.04 


8-15 


4-14 


14/16 


0.688(18) 


" 


0.05 


8-16 


4-15 


18/18 


0.696(14) 


" 


0.06 


8-16 


4-15 


18/18 


0.707(10) 


" 


0.08 


8-16 


4-15 


17/18 


0.725(7) 


" 


0.10 


8-16 


4-15 


16/18 


0.740(5) 


0.10 


0.04 


9-16 


6-15 


21/15 


0.716(9) 




0.05 


9-16 


6-15 


23/15 


0.722(7) 




0.06 


9-16 


6-15 


25/15 


0.729(6) 




0.08 


9-16 


6-15 


29/15 


0.744(4) 



TABLE III: Selected results (those used for the quark 
mass extrapolation) of the three-parameter fits for the zero- 
momentum correlators. 



fy / fy are then obtained from 



Ra(0) 



fv 



Rjip) 



PiU 



"v fv 



(15) 



For the vector mass mv in the second equation we may 
use the value obtained from our fits to the same correla- 
tors; i.e., m v = E v (p) — p 2 . Note that here f v and 
fy are bare lattice quantities which need to be renormal- 
ized. 



IV. RESULTS 

In order to extract the appropriate values for the am- 
plitudes, A vv and A VT , from the two separate correla- 
tors, we construct a x 2 measure of the form 

1, j t*,t j 

x S-^t 1 ,^) [C J (p;t J ) - C- J xp (p;t J )} 
= F{A vv {p),A v J{p),E v {p)) , (16) 

where the "expected" values, C^, xp (p; t), are given by the 
appropriate expressions in Eqs. i|13|) and S^ 1 is the in- 
verse of the covariance matrix. We extract the three pa- 
rameters by minimizing the above function, taking care 
to use large enough minimum t 1 values, thereby avoiding 
excited-state contributions. See Table ITTll for the results 
of the zero- momentum fits. 

Since we are interested in the ratios of the resulting am- 
plitudes, we repeat these fits within a single-elimination 
jackknife routine, providing the errors for fy / f v (see the 
final column of Table IIII|) . 

Figure displays the ratios of the vector meson cou- 
plings as a function of the dimensionlcss quantity, r^ni^g. 
For a wide range of quark masses the data display only 
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FIG. 2: Ratio of vector meson couplings (at the lattice scale, 
/i = 1/fl) as a function of rQrn% s for the two lattice spac- 
ings. The dashed vertical lines denote the r\rr^ and r\m\^ 
(m,gs = 2m|c — m^) physical values. The bursts denote the 
mass interpolations at m 2 si and extrapolations to m\. 



a slight curvature and are quite compatible with a lin- 
ear behavior. Only for the smallest quark masses do the 
data move upwards. Simulations on smaller physical vol- 
umes (8 3 x 24 and 12 3 x 24) at a fixed lattice spacing (a 
= 0.15 fm) suggest that the upward trend of the ratio 
at small quark masses is a finite-size effect. This inter- 
pretation is supported by the observation that also in 
spectroscopy calculations of nucleons and their excita- 
tions the finite size effects set in at the same values of 
the quark mass [2^, |2(| ^ • In order to avoid including 
such effects, we perform the (fully correlated) extrapola- 
tion to light quark mass without some of the lightest bare 
quark masses. This is necessary only for our fine lattices 
(/3 = 8.35), where the three smallest quark masses give 
rise to pseudoscalars with mpsL < 4. As can be seen in 
the plot, it is exactly these three data points which show 
finite-size effects and we therefore exclude them from the 
chiral extrapolation. For our larger lattices (/? = 7.90), 
we always have mpgL > 4; no finite-size effects are vis- 





a = 0.15 fm 


a = 0.10 fm 


(!f\ (2G eV) 


0.801(7) 


0.780(8) 


Wms 






(£)_(2 GeV) 

\ p / MS 


0.720(25) 


0.742(14) 



FIG. 3: Similar plot to that of Fig. [5] this time for the finite- 
momentum ratio m v Rj(p)/pj, where \p\ = Pj = 2ir/L. The 
symbols have the same meaning as before. 



ible and all data points are included in the fit. On the 
heavy end, all masses with r^mpg > 5 are also excluded 
from the extrapolation. To obtain the data at the strange 
quark mass we interpolate the neighboring data linearly. 

In Table IIVI we present our final results for the ra- 
tios of the couplings, matched to the MS scheme and 
evolved to the scale of 2 GeV. These results are obtained 
from the zero-momentum correlator ratios, i?4.(0) x 
Z T /Z v (2GeV). 

Since we include no improvement in our current op- 
erators, we may expect a lattice spacing dependence of 
0(a). However, since we only have results for two val- 
ues of the lattice spacing, a continuum extrapolation be- 
comes problematic. We do note that the two results for 
the p meson are consistent with each other and they are 
also in agreement with the results of other lattice calcu- 
lations [HHi: 0.72(2)(±^) in Ref. [IJ. The trend of the 
<f> results suggests a continuum value below our (3 = 8.35 
result of 0.78, in rough agreement with the 0.76(1) result 
of Ref. |3. 





a = 0.15 fm 


a = 0.10 fm 


(!t\ (2 GeV) 

Wms 


0.805(15) 


0.746(24) 


(£)_(2 GeV) 

\ P J MS 


0.740(39) 


0.700(38) 



TABLE IV: fv/fv values in the MS scheme at p, = 2 GeV TABLE V: fv/fv values in the MS scheme at u = 2 GeV as 
as determined from the lattice ratio i?4(0). determined from the lattice ratio Rj(2n / L). 
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If we normalize using the experimental value of f°± p ' ~ 
208 MeV ITol- I20I w e findgpod agreement with the QCD 
sum rule result [lJlllEi: fplfp^' = 0.74(3). Such 
comparisons, however, are difficult to assess due to the 
different systematics in the two approaches (e.g., our lat- 
tice calculation is quenched). 

The results for the finite-momentum correlator ratio, 
Rj(p), are shown in Fig. [3] The quark-mass interpola- 
tions and extrapolations are performed just as before for 
the zero-momentum ratios. 

Tabic displays the renormalized coupling ratios ob- 
tained from the non-zero-momentum correlator ratios. 
The errors are significantly larger for these than those 
for the zero-momentum results. However, the results of 
this consistency check are compatible with those from the 
zero-momentum ratios (< 1.5ct even for the jackknifcd 
difference at and below the strange-quark mass). 

Let us briefly summarize our findings. We calculate 
fy/fy going considerably closer to the chiral limit (at 
least on our larger lattices where the results do not dis- 
play significant finite-volume effects) than previous cal- 
culations. This allows us to monitor the mass depen- 



dence of this ratio of couplings over a much larger mass 
range. We find the quark mass dependence to be rela- 
tively weak, as in QCD sum rule calculations. The val- 
ues we obtain agree well with the extrapolated values 
of other lattice calculations at larger quark masses. We 
include a consistency check using spatial tensor corre- 
lators at finite momentum and find compatible results. 
Our final numbers are f£/f p = 0.720(25), 0.742(14) and 
/J-// = 0.801(7), 0.780(8) at a = 0.15, 0.10 fm, respec- 
tively. 
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